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Long-term ambient concentrations of inhalable particles less than 10 pm in diameter (PM 10 ) (1973- 
1992) and other air pollutants—total suspended sulfates, sulfur dioxide, ozone (0 2 ), and nitrogen di¬ 
oxide—were related to 1977-1992 mortality in a cohort of 6,338 nonsmoking California Seventh-day 
Adventists. In both sexes, PM 10 showed a strong association with mortality for any mention of nonma- 
lignant respiratory disease on the death certificate, adjusting for a wide range of potentially con¬ 
founding factors, including occupational and indoor sources of air pollutants. The adjusted relative 
risk (RR) for this cause of death as associated with an interquartile range (IQR) difference of 43 d/yr 
when PM 10 exceeded 100 pg/m 3 was 1.18 (95% confidence interval [Cl]; 1.02, 1.36). In males, PM 10 
showed a strong association with lung cancer deaths—RR for an IQR was 2.38 (95% Cl: 1.42, 3.97). 
Ozone showed an even stronger association with lung cancer mortality for males with an RR of 4.19 
(95% Cl: 1.81, 9.69) for the IQR difference of 551 h/yr when 0 3 exceeded 100 parts per billion. Sulfur 
dioxide showed strong associations with lung cancer mortality for both sexes. Other pollutants 
showed weak or no association with mortality. Abbey DE, Nishino N, McDonnell WF, Burchette 
HI, Knutsen SF, Beeson WL, Yang IX. Long-term inhalable particles and other air pollutants 
related to mortality in nonsmokers. am jrespir crit care mej> isssn53^73-3*2. 


There have been three recent prospective cohort studies in the 
United States of the relationship between long-term ambient 
concentrations of particulate air pollution and mortality (1-3). 
The Six Cities Study followed a cohort of 8.111 adults, ran¬ 
domly sampled between 1974 and 1977 from six cities located 
in the eastern and mid western United States (1). The Ameri¬ 
can Cancer Society (ACS) Study followed a cohort of 552,138 
adults who were enrolled by ACS volunteers throughout the 
United States and Puerto Rico in 1982 (2) A third study, 
known as the Adventist Health Study of Smog (AHSMOG 
Study), has followed a cohort of 6,338 nonsmoking California 
Seventh-day Adventists (SDAs) since 1977 (3). Two of these 
studies—the Six Cities Study and the ACS Study—have shown 
positive associations between long-term concentrations of 
ambient particles and increased risks of mortality (1, 2). The 
AHSMOG Study, however, showed no such associations us¬ 
ing mortality data from 1977-1987 (3). 

Mortality ascertainment has recently been updated on the 
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AHSMOG cohort through 1992. Previously, only indirect esti¬ 
mates of inhalable particles less than 10 pm in diameter 
(PM W ) formed from site- and season-specific regression equa¬ 
tions on total suspended particles (TSP) were available (4). 
Since 1987, PM 1C > has been monitored throughout California. 
These data have been used to update ambient concentration 
estimates for study participants using actual monitored PM 10 
through 1992. The purpose of this paper Is to study associa¬ 
tions between mortality and long-term ambient concentra¬ 
tions of PM,, and other air pollutants using the recently up¬ 
dated mortality data through 1992. 

METHODS 

Questionnaire Data and Follow-up of the Cohort 

The study methods have been described in detail previously (5-7). 
Briefly, 6,338 nonsmoking. non-Hispanic white SDA residents of Cal¬ 
ifornia were enrolled in 1977. The study participants, ages 27-95 yr at 
baseline, were part of a larger Adventist Health Study (AHS) (8). in¬ 
clusion criteria were the following: having lived 10 yr or longer within 
5 miles of their residence at time of enrollment; and residing in one of 
the three California air basins of San Francisco, South Coast (Los An¬ 
geles and eastward), or San Diego, or being part ol a 10% random 
sample of AHS subjects from the remainder of California who met 
the other inclusion criteria. 

All subjects completed a detailed lifestyle questionnaire as part of 
the AHS in 1976. This questionnaire ascertained anthropometric 
{height, weight) data, exercise patterns, use of alcohol and tobacco, 
occupation, current and past dietary habits, parental history of cancer, 
and history of selected medical conditions (8). Monthly residence artd 
work location histories were obtained for each study subject for the 
period 1960 through April 1992 or date of death (if subject died prior 
lo April I, 1992) by using mailed questionnaires 11977, 1987, and 
1992), tracing by telephone, and interviewing of surrogates for de 
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ceased subjects. Only 156 (2.5%) subjects were lost to follow-up. 
These latter individuals were censored at dace of last contact for pur¬ 
poses of inclusion in risk sets. The 1977,1987, and 1992 questionnaires 
□Iso contained standardized respiratory symptoms questions, now in¬ 
cluded as part of the American Thoracic Society questionnaire (9). 
Computer algorithms were used to classify individuals as to whether 
they had definite, possible, or no symptoms of asthma, chronic bron¬ 
chitis, or emphysema (5). The questionnaires also ascertained lifestyle 
and housing characteristics pertinent to relative exposure to ambient 
air pollutants as well as indoor sources of air pollutants. 

Air Pollution Estimates 

Estimates of monthly ambient concentrations of PMjo, ozone (Oj), 
sulfur dioxide (SO.). and nitrogen dioxide (N0 2 ) were formed for 
study participants for the period 1966-1992 using fixed site monitor¬ 
ing stations maintained by the California Air Resources Board. For 
suspended sulfates (SQj). the time period was 1977-1992 since S0 4 
was not monitored in California prior to 1977. The detailed methods 
for estimating ambient air pollutants for study participants are de¬ 
scribed elsewhere (4, 6, 10), Briefly, monthly indices of ambient air 
pollutant concentrations at 348 monitoring stations throughout Cali¬ 
fornia were interpolated to zip code centroids according to home or 
work location histories of study participants, cumulated, and then av¬ 
eraged over rime. Interpolations were restricted to Zip code centroids 
within 31.25 miles (50 km) of a monitoring station and were not al¬ 
lowed to cross barriers to air flow or other topographic obstructions in 
excess of 250 m above the surrounding terrain (6). Quality ratings that 
reflected the distance of the nearest monitoring station from the zip 
code centroid were attached to all interpolations in order to allow sen¬ 
sitivity analyses (4,6). 

Concentrations of PM 1C prior to 1987 were estimated using site- 
and season-specific regressions based on TSP (4). Since 1987, PM ]C 
has replaced TSP as the particulate pollutant monitored throughout 
California, so monitored F VI ] n data were available for the later years 
of the study. For PMm and 0 3 . exceedance frequencies and excess 
concentrations above several cutoffs were estimated in addition to 
mean concentrations. Exceedance frequencies were defined as the 
sum of days above a specified cutoff of PM 10 and the sum of hours 
above a cutoff for 0 3 . The cutoffs used for PM in were 40, 50. SO. 80, 
and 100 pg/m 3 , denoted by PM lo (40), etc. The cutoffs used For 0 3 
were 60, 80,100,120, and 150 parts perbillion (ppb). An additional in¬ 
dex used for 0 3 was the monthly average of the daily 8-h average from 
9:00 a.m. to 5:00 p.m. (used to correspond to the usual hours at work 
locations because separate interpolations were used for work loca¬ 
tions). Excess concentration indices were also computed with excess 
concentration defined as the integrated sum of time X concentration 
in excess of a cutoff. Correlations of these indices with the exceedance 
frequency indices were close to one for corresponding cutoffs, so the 
results are. not reported separately for them (6, 11). 

Mortality Ascertainment 

Deaths during the follow-up period were ascertained using three 
methods: (7) computer-assisted record linkage wiLh the California 
death certificate files for the years 1977-1992; (.5) National Death In¬ 
dex files, 1979-1992; and {3; our tracing procedures, which included 
church records. We thus identified 1,628 (989 females, 639 males) 
study subjects who died. Death certificates were obtained for all of 
these, and a state-certified nosalogist blinded to exposure data coded 
each death certificate according to the International Classification of 
Disease. Ninth Revision (1CD9) codes 

Statistical Analysis 

Sex-specific adjusted mortality relative risks (RRs) were estimated us¬ 
ing Cox proportional hazards regression models. Combined male and 
female RRs were also estimated, stratifying on sex. when sex-specific 
models had similar regression coefficients for covariates. Time-depen- 
dent proportional hazards models were used that averaged the pollut¬ 
ant from the first monthly vaiue through the month of death. For lung 
cancer mortality, a 3-yr lag was used, which averaged the pollutant 
only up to the full 3 yr prior lo the month of death because of the ex¬ 
pected long latency period between exposure and incidence of cancer. 
This is the maximum lag time that can be used with the present data. 


since monthly values prior to 1973 were not stored in the data. The 
time variable used in the models was survival time from date of enroll 
ment, except for lung cancer mortality, where age on study was used 
due to the expected lack of short-term effects. Using time on study 
rather than age as the time variable can bias proportional hazards re¬ 
gression models unless none of the covariates in the models are corre¬ 
lated with age at baseline (12)- We found this condition to be met for 
models in which we used time on study with all correlations having an 
ahsolute value less than 0.17. Survival times of participants who did 
not die were censored at the end of the study period or at the day of 
last contact if the participant was lost to follow-up. Adjusted RRs 
were calculated and reported for differences in air pollutant concen¬ 
trations equal to the interquartile ranges (IQRs) of the 1973-1992 av¬ 
erages of pollutants observed across the study population. Proportional 
hazard models were stratified if needed to satisfy the proportional haz¬ 
ards assumption. Failure of the proportional hazards assumption was 
indicated if the statistical significance for change in log likelihood of 
the fit of the model was less than 0.05 for the cross-product of any 
variable in the model with the log of the time variable (13). 

The following covariates were controlled for in a priori models: 
age at baseline. April 1, 1977; pack-years of past smoking; years lived 
with a smoker; years worked with a smoker; years of education at 
baseline; an indicator variable for employment in an occupation with 
significant exposure to air pollutants for more than 10 yr (14); and 
body mass index (BMI) measured as {weight in kg)/(height in nil" 
Body mass index was coded as a series of indicator variables to allow 
for nonlinear relationships with mortality. For lung cancer mortality, 
the only covariates used were age, pack-years of past smoking, and 
education due to the small number of cases. Separate Cox propor¬ 
tional hazards models were estimated for the following cause of death 
categories; all natural cause (ICD9:001-799), cardiopulmonary dis¬ 
ease (ICD9:40l-44Q and 460-519). nonmalignant respiratory disease 
(ICD9 460-519), and lung cancer (1CD9T62). Only the underlying 
cause of death was used. An additional outcome used for nonmalig¬ 
nant respiratory disease was any mention of nonmalignant respiratory 
disease as an underlying or contributing cause of death, subsequently 
abbreviated CRC mortality. Deaths not in a specific category were 
censored at time of death. Robustness of the models was evaluated by 
estimating each model separately for never smokers versus past smok¬ 
ers as well as conducting subgroup analyses for other key variables. 

In order to further ride out confounding or effect modification for 
variables not mentioned above, additional models were developed for 
those causes of death that showed strong positive associations with 
pollutants—all natural cause, any mention of nonmalignant respira¬ 
tory disease, and lung cancer. A large number of potentially con 
founding Or effect-modifying factors. Such as medical history, dietary, 
or other lifestyle factors, were evaluated to be sure their inclusion did 
not substantially change the pollutant effect estimate. The Appendix 
gives a detailed description of these factors. Those factors were in¬ 
cluded that substantially changed the pollutant regression coefficient 
when included or that significantly (p < 0.05) improved the fit of the 
model according to the likelihood ratio criteria (15). A change in the 
PMmQOO) regression coefficient was defined as substantial if it re¬ 
sulted in a change in the RR corresponding to an IQR difference of 
PM lo (100) oTmore than 0.D4. To assure parsimonious models, covari¬ 
ates that did not meet these criteria were dropped. These refined 
models were used for most results and sensitivity analyses. Regression 
diagnostics were conducted for all refined models. These regression 
diagnostics included checks of the proportioned hazards assumption, 
checks of log linearity for continuous variables, checks for outliers, 
and stratified analyses. Two-pollutant models were used to evaluate 
the stability of single pollutant associations. 

RESULTS 

Table 1 shows the numbers of sex-specific deaths occurring In 
the cohort between 1977 and 1992. Descriptive characteristics 
of the cohort at baseline in 1977 are included for both sexes. 
Table 2 gives descriptive statistics for ambient air pollutants 
estimated for study participants averaged 1973 through month 
of censoring Table 3 shows adjusted RRs for different causes 
of mortality. The RRs are calculated for a difference in the 
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TABLE 1 

NUMBER Of CAUSE-SPECIFIC DEATHS AND BASELINE (1977) CHARACTERISTICS 
FOR COHORT (4,069 FEMALES, 2,278 MALES) 


Female 

Male 


989 

639 

Total deaths 1977-1992 

965 

610 

All natural Cause deaths (ICD9:0D1-799) 

631 

398 

Cardiopulmonary deaths (ICD9':401-440and 460-519) 

72 

63 

Nonmalignant respiratory deaths (1CD9:460-519) 

246 

164 

Deaths with any mention of nonmalignant respiratory disease in underlying or contributing cause of death* 

12 

18 

Lung cancer deaths (ICD9:162) 

59.2 

58.5 

Age in years, mean 

13.8 

36.3 

% Smoked in the past 

11.3 

19.8 

Pack-years of smoking, mean for past smokers 

47.8 

33.8 

% Ever lived with smoker 

19.9 

16.1 

Years lived with smoker, meant 

38.2 

48.0 

% Ever worked with smoker 

11.3 

15.6 

Years worked with smoker, mean' 

0.4 

10,3 

% Occupational exposure to air pollutants for more than lOyr 

13.1 

14,3 

Years of education. mean 

24.6 

25 0 

Body mass index, (weight in kg)/fheight in m) 2 , mean 

6.4 

9.5 

% Currently use alcoholic beverages 

35.2 

44.3 

% High total exercise level 

9.0 

17.5 

Hours outdoors per week, mean 

33.9 

298 

% Whose father or mother had cancer 

32.3 

27.6 

% High fruit/vegetabte consumption' 

46.6 

40.0 

% High antioxidant vitamin consumption From pills 5 

33.0 

28.8 

% Prior heart attack, stroke, diabetes, or high blood pressure 


* All death counts except this one are for underlying cause only 

* Mean of non-zero values. 

1 Two or more servings of fruit or fruit juice p-er clay and two or more servings of cooked green vegetables, green salads, or tomatoes pier 
day, 

s At least daily consumption of vitamin A, or more than 1,000 mg vitamin C per week, or 200 U or more vitamin E per week, or 1 or 2 
multivitamin pills daily 


IQR of the pollutant across the study population. Table 4 
shows RRs for the other covariates in sex-specific refined 
models for CRC and lung cancer mortality. These were the 
mortality outcomes showing strong associations with air pol¬ 
lutants. Refined and a priori models gave similar results for 
these outcomes. 

All Natural Cause Mortality 

For males, PM 10 showed a significant association with all natu¬ 
ral cause mortality, but none of the other pollutants did. When 
other pollutants were added one at a time to form two-pollut- 
ant models with PM I0 , the PM^ coefficient remained stable. 
When any mention of malignant or nonmalignant respiratory 
mortality was removed from all natural cause mortality, the 


RR associated with the IQR difference in PM ]o (10O) dropped 
from 1.12 to 1.09 (95% confidence interval [Cl]: 0.96, 1.23). 

For females, none of the pollutants showed a positive asso¬ 
ciation with all natural cause mortality. 

Cardiopulmonary Mortality as Underlying Cause 

None of the pollutants showed a significant association with 

cardiopulmonary' mortality for either sex. 

Nonmalignant Respiratory Mortality as 
Underlying Cause 

None of the pollutants showed a significant association with 
nonmalignant respiratory mortality for either sex. It was noted 
that the proportional hazards assumption was not satisfied 


TABLE 2 

DESCRIPTIVE STATISTICS TOR ESTIMATED AVERAGE AMBIENT AIR POLLUTANT 


VALUES 1973 TO 

CENSOR 

DATE' FOR 

SUBJECTS 



Pollutant 

n* 

Mean 

Standard 

Deviation 

Range* 

Interquartile 

Range 

PM 10 mean concentration in p.g/m 3 

5,963 

51.24 

16.63 

B3.93 

24.08 

Days/year above 100 p-g/m 3 PM 10 

5,991 

31.08 

32.48 

173.84 

42.63 

Suspended sulfate (S0 4 ) mean concentration in i-tg/m 3 

5.070 

7.24 

2.55 

32.11 

2.97 

Sulfur dioxide (S0 2 ) mean concentration in ppb 

4.353 

5.62 

2.81 

18.96 

3.72 

Ozone (0 3 ) mean concentration in ppb 

5,893 

26.11 

7.65 

43.91 

12 03 

Hours/year above 100 ppb 0 3 

5,893 

329.61 

294.51 

987,97 

551.10 

Nitrogen dioxide (NO z ) in ppb 

5,652 

36.76 

12.99 

67.67 

19.78 


Definition of abbreviation'. PfVl,^ = inhalable particufatErs less than 10 p.m. 

* Except S0 4 , 197 7 to censor date. 
f Number of subjects with 80% non-missing data. 

1 Maximum average value — minimum average value with average computed over entire time period. All minimum values were 2 ern so 
that the ranges presented are maximum values, The zero values v/ere assigned to a few individuals living in pristine areas remote from 
monitoring stations. 
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TABLE 3 

AD1USTED MORTALITY RELATIVE RISKS AND 95% CONFIDENCE INTERVALS (IN PARENTHESES) BY 
CAUSE OF DEATH FOR AN INTERQUARTILE RANGE {IQR) DIFFERENCE IN SPECIFIC AIR POLLUTANTS 





Pollutant Average 1973-Censor Date 




PM,o above 

PM, 3 Mean 

S0 4 Mean 

S0 2 Mean 

0 3 above 

0 3 Mean 

N0 2 Mean 


TOO |ig/m 3 

Cone 

Cone.' 

Cone. 

100 ppb 

Cone. 

Cone. 

IQR 

Cause of death 

$3 d/yr 

24.08 |xg/m 3 

2.97 pg/rrr 

3.72 ppb 

551.1 h/yr 

12.03 ppb 

19.78 ppb 

All natural cause 








Males* 

1.12** 

1.11 

1.05 

1.05 

1.14 

1.09 

1.03 


(1.01, 1.24) 

(0.98, 1.26) 

(0.95, 1.16) 

(0.94, 1.18) 

(0.98, 1.32) 

(0.95. 1.25) 

(0.91, 1.17) 

Females* 

094 

0.94 

0.94 

1.00 

0.90 

0.95 

0.99 


(0.86, 1.03) 

(0.84, 1.04) 

(0.86, 1.02) 

(0.91, 1.10) 

(0.80, 1.02) 

(0.85. 1 06) 

(0.89, 1.11) 

Cardiopulmonary 








Males* 

1.09 

1.10 

1-05 

1.01 

1.06 

1.O8 

1.01 

F&males* 

(0 95. 1.24) 

(0.94, 1.30) 

(0.28. 3.94) 

(0.87, 1.18) 

(0.88, 1.29) 

(0.91.1.29) 

(0.86, 1.19) 

0 90 

0.92 

0.97 

1.02 

088 

0.97 

1.03 


(0.80. 1.01) 

(0.80,1.05) 

(0.87. 1.08) 

(0.90. 1.15) 

(0.75. 1.02) 

(0.84. 1.12) 

(D.90, 1.18) 

Any mention of nonmalignant 








respiratory 

Males 6 

1.28*’ 

1.23 

1.04 

087 

1.20 

1.12 

0.87 


(1,03, 1,57) 

(0.94, 1.61) 

(0.84, 1,29) 

(0.68, 1,11) 

(0.88, 1.64) 

(0 85. 1.47) 

(0.67, 1.14) 

Females* 

1.10 

1.10 

1.02 

0.98 

1.01 

1.05 

0.97 


(0.91 r 1.33) 

(0.86 r 1.40) 

(0.84, 1.25) 

(0.79. 1.23) 

(0.77, 1.33) 

(0.82, 1.35) 

(0.76, 1.23) 

Lung cancer 








Males 1 

2.38" 

3.36 Tr 

— 

1.99 ,r 

4 19” 

2.10 

1.82 


(1.42, 3.97) 

(1.57,7.19) 


(1.24. 3.20) 

(1.81. 9.69) 

(0.99, 4.44) 

(0.93, 3.57) 

Females* 

1,08 

1.33 

— 

3.01" 

1.39 

0.77 

2.81** 


(0.55. 2.13) 

(0.60, 2.96) 


(1.88. 4.84) 

(0.52, 3.67) 

(0.37, 1.61) 

(1.15, 6.89) 


Definition of abbreviations: cone. - concentrations; PM 10 = inhalable particulates less than 10 jun. 

4 Averaged 1977-censor date. 

1 Covarialea included age, years of education, pack-years of past smoking, history of high blood pressure, years lived with a smoker, and 
total exercise level. 

* Covariates included age,, years of education, pack-years of past smoking, years lived with a smoker, years worked with a smoker, ocraj- 
pational exposure for more than 10 yr, and body mass index as two dummy variables with reference category between stated values (< 21,0 
or ^ 27 5 for females. < 22.4 or ^ 27.3 for males). 

* Covariates included patik-years of past smoking, body mass index as two dummy variables with reference category between stated val¬ 
ues (< 21.0 cr > 27-5 for females, < 22.4 or 27.3 for males), total physical exercise, and age within age strata Models were stratified by 
age (< 65 versus > 65) andl excluded over age 85 in 1977. 

1 Covariaccs included years of education, pack-years of past smoking, and alcohol use, 

r Covariates included yealrs Of education and pack-years of past smoking. 

** p < 0.05. 

"pc 0.01. 

•*p < 0.001. 


with the a priori model. A reduced model that did satisfy the 
proportional hazards assumption found a positive but nonsig¬ 
nificant association with PM| 0 (100) for males, whose RR for 
an IQR Increase in PMic(lOO) was 1.15 (95% Cl: 0.81, 1.63). 

Any Mention of Nonmalignant Respiratory Mortality 

Only PM m showed a significant association with any mention 
of nonmalignant respiratory mortality. Since sex-specific mod 
els for this outcome gave similar regression coefficients for 
PM 1U and other covariates (Table 4), combined male and fe¬ 
male RRs for PM|„ were also calculated. The combined male 
and female adjusted RRs of CRC mortality for the IQRs of 
PMi 0 (100) and PM 10 mean concentration, respectively, were 
1.18 (95% Cl: 1.02, 1.36) and 1.16 (95% Cl: 0.97, 1.39). Esti¬ 
mation of the RRs for CRC mortality used stratification by 
age and sex in order to satisfy the Cox proportional hazards 
assumption. It was necessary to exclude subjects over 85 yr of 
age in 1977 from the proportional hazards regression models 
since the proportional hazards assumption could not be satis¬ 
fied when this older age group was included No association 
was observed between CRC mortality and PMio for the 151 
male and female subjects over age 85 In 1977 When other pol¬ 
lutants were added one at a time to form two-pollutant models 
with PMjfl, the P\l ]0 coefficient remained stable or increased. 


Lung Cancer Mortality (Males) 

In single-pollutant analyses, PM m (100), Oj(1 00), and S0 2 
were all significantly associated with increased risk of lung 
cancer mortality in males (Table 3). When N0 2 was added 
with these pollutants in turn to form two-pollutant models, 
their coefficients remained stable and N'O z remained nonsig¬ 
nificant. When S0 2 was added to the PIvl 10 and 0 3 models, all 
coefficients remained stable, suggesting that the effects of S0 2 
are independent of these other pollutants, as is also suggested 
by the low correlations of SO z with other pollutants (Table 5). 
It was difficult to separate the effects of PMio and 0 3 because 
they were much more highly correlated. In two-pollutant 
models of PM ]o (10O) and O 3 (100), both coefficients remained 
positive and were reduced in magnitude, with the Qj coeffi¬ 
cient remaining stronger than the PM,„ coefficient. 

Lung Cancer Mortality (Females) 

Both SO; and NO; were significantly associated with in¬ 
creased risk of lung cancer mortality in females (Table 3); 
however, in a two-pollutant mortality model with S0 3 , the 
N0 2 coefficient became negative, while the SO, coefficient re¬ 
mained strongly positive. The SO z coefficient also remained 
stable in two-pollutant models with PM 1lf (100) and 0 3 (100). 
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TABLE 4 


SEX-SPECIFIC ADJUSTED MORTALITY RELATIVE RISKS AMD 95% CONFIDENCE INTERVALS 
(IN PARENTHESES) ASSOCIATED WITH INCREMENTS OF PM 10 C100) AND OTHER COVARIATES 
FOR SELECTED MORTALITY OUTCOMES (4,060 FEMALES, 3,778 MALES) 



CRC Mortality* 

Lung Car>cer Mortality 

Covariate 

Females 

Mates 

Females 

Males 

Number of cases 

152* 

1 20* 

12 

17 

43 d/yr when PM m exceeded 100 ug/m 5 

1.10 

1.28" 

1.09 

2.38' 


(0.91, 1.33) 

(1.03, 1.57) 

(0 55. 2.13) 

{1.42, 3.97) 

5 yr age 

i.ao’ 

(1.56, 2.08) 

2,12’ 

{1.79. 2.52) 



4 yr of education 

— 

— 

0,89 

(0.39, 2.00) 

0.58 

(0.30, 1.11) 

10 pack-years of past smoking r 

1,14 

1.20* 

1.41 

1.39* 


(0.86, 1.49) 

(1.12. 1.29) 

(0.77. 2.57) 

(1.19, 1.62) 

1 unit increase in total exercise level* 

0.86 1 

(0.76, 0.98) 

0.91 

(0.77, 1.09) 

— 

— 

Lowest BMI category* {females < 21.0, males < 22.4) 

1.98' 

(1.34. 2 94) 

1.20 

(0.76,1.83) 

— 

— 

Highest BMI category 4 [females & 27.S, males > 27.3) 

2,11* 

[1,44, 3.11) 

1.10 

(0.66, 1.83) 

— 

— 

Current alcohol use in 1977 

— 

— 

— 

3.g2i' 

(1 35. 11.43) 


Definition of abbreviations: BMI - body mass Endex: CRC = any mention of nonmallgnanl respiratory mortality ir underlying or contrib¬ 
uting cause of death; PM 10 - irihafable particulates less than 10 p-m. 

* Stratified by age. excluded over age 95 in 1977, age strata were < 65 vs. s=> 65. 

' Mean years since last smoked: 24.5 for females, 28.4 for mates. 

1 Four levels: 1 = light/none; 2 = low; 3 = moderate; 4 = high. 

5 Reference categories: > 21.0 and < 27.5 for females; ^ 22.4 and < 27.3 for males. 

^ p < 0.05. 

' p < 0.001. 

M Age used as the time variable. 


Robustness or Associations for Subgroups 

We investigated the robustness of the observed air pollutant 
associations by conducting subgroup analyses for CRC mor¬ 
tality and lung cancer mortality for the pollutant that had the 
strongest association with an outcome—PM ]o {100) for CRC 
mortality, O a for lung cancer mortality in males, and S0 2 for 
lung cancer mortality in females. The combined sex model was 
used for CRC mortality, 

The association of CRC mortality with PM 10 was similar in 
past smokers. RR = 1.23 (95% Cl: 0.90, 1.66), and never 
smokers, RR = 1.17 (95% Cl; 1.00, 1.37) for an 1QR increase 
of PM| O (100). The associations of lung cancer mortality with 
0 :l in males were strongly positive for both past smokers. 
RR — 4.25 (95% Cl; 1.50, 12.07), and never smokers, RR = 
6.94 (95% Cl: 1.12, 43.08), for the IQR of O s . For female 


TABLE 5 

CORRELATIONS BETWEEN LONG-TERM AVERAGES OF POLLUTANTS 
ESTIMATED FOR STUDY PARTICIPANTS, 1973 THROUGH MONTH 
OF CENSORING (EXCEPT SO„, ONLY SINCE 1977) 



PM 10 me 

Ozone (100) 

Ozone me 

SO* me 

S0 2 me 

N0 2 me 

PM 1D <100) 

0.85 

0.84 

0.63 

0.33 

-0.05 

0.15 

PM 10 me 

1.00 

0.83 

0.77 

0.68 

0.31 

0.56 

Ozone 100 

— 

1.00 

0.78 

0.43 

0.13 

0.41 

Ozone me 

— 

— 

1.00 

0.53 

0.09 

0.36 

SO, me 

— 

— 

— 

1.00 

0.68 

0.76 

SO, me 

— 

- 

- 

- 

1.00 

0.79 


Definition of abbreviations'. PM n0 - inrialabte particulates less than 10 urn; me - 
mean concentration: PMi O (1O0) = d/yr when PM 10 exceeded 100 pig/m 3 ; ozone (100) = 
h/yr when ozone exceeded 100 ppb 

* Number of subjects having al least 80% noTimissing monthly values tor both pol¬ 
lutants varies from 4,150 to 5,963. 


never smokers, an RR of 2.99 (95% Cl: 1.66, 5.40) was ob¬ 
served for the IQR of SOj. There were too few female past 
smokers to calculate the RR. We did not have a direct mea¬ 
sure of years since last smoked for past smokers. However, 
years since joining the church could be taken as a crude lower 
bound on this variable. Of past smokers, 7.3% had joined the 
church less than 5 yr ago and could thus be classified as possi¬ 
ble recent past smokers. When these individuals were ex¬ 
cluded from analyses, the regression coefficient of PM 10 for 
CRC mortality in both males and females and that of O3 for 
lung cancer mortality in males became a little larger, and the 
coefficient for lung cancer mortality and SO, In females, 10% 
smaller. When dummy variables representing years since join¬ 
ing the church (0-16 yr as a reference group, 17-33 yr, and 34- 
76 yr) were added to the model of CRC mortality and PMj O (100) 
for past smokers, the regression coefficients for PM 10 and smoke 
pack-years were within 12% of previous values. 

When individuals who reported current use of alcoholic 
beverages at baseline (6.4% of females and 9.5% of males) were 
excluded, stronger associations were observed for PMjo(lOC) 
with CRC mortality, 0 3 (100) with lung cancer mortality in 
males, and S0 2 with lung cancer mortality in females. 

For 2,680 subjects who reported low antioxidant vitamin 
consumption at baseline, the RR for CRC mortality for an 
IQR difference of PM I0 was 1.26 (95% Cl: 1.01, 1.57), while 
the RR for 2,326 subjects with high antioxidant vitamin con¬ 
sumption was 1.08 (95% Cl: 0.89, 1.32). For lung cancer mor¬ 
tality, there were too few cases to conduct subgroup analyses. 

For subjects who reported a prior heart attack, stroke, dia¬ 
betes, or high blood pressure at baseline, the association of 
PM I() with CRC mortality was similar to those who did not 
report these events. In 1977, subjects were classified as to 
whether they had definite symptoms of asthma, chronic bron- 
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chitis, or emphysema (5). The presence of these respiratory 
symptoms in 1977 increased the risk of death from respira¬ 
tory-associated causes to 1.49 (95% Cl: 1.12, 1.99). Subjects 
classified as having definite respiratory symptoms in 1977 
showed a similar association of PM 1(] with CRC mortality to 
those without such symptoms. 

Subjects were stratified into three groups according to 
hours per week spent outdoors as reported at baseline: low 
was classified as 4 h/wk or less, medium as greater than 4 but 
less than or equal to 16 h/wk, and high as more than 16 h/wk. 
The RRs of CRC mortality for the IQR of PM| 0 (100) for each 
of these groups were; low, RR = 1 07 (95% Cl; 0,85, 1.34); 
medium, RR = 1.18 (95% Cl; 0.90,1.55); and high, RR = 1.31 
(95% Cl: 1.02, I 68) 

When subjects were restricted to those living in high popu¬ 
lation density areas (> 10 homes within a quarter-mile ra¬ 
dius) , results were similar to those reported for PM IC , and CRC 
mortality in both sexes. For lung cancer mortality in males, the 
RR as associated with elevated 0 3 was higher. This was true 
for PM 10 as well. 

Potentially Confounding Factors 

A large number of lifestyle and other potentially confounding 
factors of PM (0 were evaluated in this study (Table 1 and Ap¬ 
pendix). None of these potentially confounding factors, when 
added one at a time to the PM] 9 mortality models, substan¬ 
tially changed the PM 10 coefficient, except possibly for factors 
already included in refined models; thus, adjustments for these 
factors were not needed. We defined a substantial change in the 
PM ]a (100) regression coefficient as one that would result in a 
change of more than 0.04 in the RR corresponding to an IQR 
difference of PM lf) (100). The effects on the PM,,, coefficient of 
the potential confounders that were included in the refined 
sex combined CRC mortality model were as follows. For a 
basic model stratified by sex and age and including only age as 
a covariate, the RR of CRC mortality for an IQR increase in 
PM lo (100) was 1.10 (95% Cl: 0 97, 1.25); when smoke pack- 
years were added, the RR was 1 14 (95% Cl: 1.00, 1.30); 
when the BMI indicator variables were also added, the RR 
was 1.18 (95% Cl: 1.02, 1.35); with total exercise scale added 
to the above variables in the model, the RR was 1.18 (95% Cl: 
1.02, 1.36). 

To further evaluate potential confounding, we conducted 
stratified analysis for any potential confounders that showed a 
strong association with PM 10 and a stronger positive associa¬ 
tion with CRC mortality than PM l0 . An association of a poten¬ 
tial confounder with PM ln was deemed strong” for a continu¬ 
ous variable if the absolute value of the correlation coefficient 
with PM lo (10O) was greater than 0.20. For categorical vari¬ 
ables, an association was deemed "strong" if the p value of a 
one-way analysis of variance of PM, rl (U)0) across the catego¬ 
ries was less than 0 10. An association of a potential con- 
founder with CRC mortality was deemed stronger than that 
for PM| 0 (100) if the magnitude of the RR for the potential con- 
founder exceeded the RR for the IQR increase of PM 10 (100). 
For continuous valued potential confounders. we evaluated 
the RR for an IQR increase: for categorical potential con¬ 
founders, we evaluated the RR with reference to the lowest or 
highest category. Based on these criteria, stratified analyses 
were conducted for both sexes; occupational exposures to air 
pollutants for more than 10 yr, current alcohol consumption 
(yes, no), vitamin C intake, multivitamin intake, a salad index, 
and several meat indices, in the slratified analyses, we evalu¬ 
ated the RR of CRC mortality as associated with the IQR in¬ 
crease of PiVf ut (100) within each stratum of the confounding 
factor. The RRs were consistent with those reported above for 


the entire group of subjects, except for some small subgroups, 
suggesting that no confounding was present. For each factor 
evaluated, CIs for all subgroups overlapped each other. A simi¬ 
lar evaluation of potential confounders was done for 0 3 and 
lung cancer mortality for males. Again no evidence of confound¬ 
ing was found. Although we evaluated a large number of po¬ 
tential confounders, we cannot rule out the possibility of con¬ 
founding by unmeasured variables. 

Possible Underreporting of Current Smoking 

All reported current smokers in 1977 were excluded from all 
analyses. However, since smoking is a religious proscription of 
the SDA Church, it is possible that current smoking was un¬ 
derreported This could bias our observed associations be¬ 
tween long-term ambient concentrations of air pollutants and 
respiratory-related mortality (16). Using the reported current 
smoking rates in our cohort in 1977 (1.1% for females, 2.2% 
for males), we conducted sensitivity analyses, assuming a 10- 
fold greater risk for current smokers and 50% underreporting, 
which was (/) nondifferential and then {2) differential with re¬ 
spect to air pollution exposure. For case (2), we assumed that 
underreporting occurred in the top quartlle of air pollution ex¬ 
posure, but not in the other quartiies. Our sensitivity analyses 
indicated that the RRs reported in Table 3 would underesti¬ 
mate the true RR for the case of nondifferential underreport¬ 
ing and would overestimate the true RRs by not more than 
15% of the true RR for the case of differential underreporting. 

Cutoffs 

The relationship between CRC mortality in both sexes and 
PMiu was evaluated for other cutoffs—40, 50, 60, and 80 p.g/ 
m a . Exceedance frequencies for cutoffs of 60 |xg/m 3 and higher 
showed statistically significant associations with CRC mortal¬ 
ity. A similar evaluation for lung cancer mortality in males and 
the different cutoffs of 0 3 showed all cutoffs of 60-150 ppb as 
well as the 8-h average to be significantly associated with lung 
cancer mortality, but not the 24-h average of 0 3 . 

Time Trends 

Except for 0 3 mean concentration, levels of all pollutants and 
pollutant indices used in this study declined from 1973 to 1992, 
or 1977 to 1992 in the case of SO., Ozone mean concentration 
stayed about the same, though exceedances above 100 ppb Oj 
declined. The most marked decline was seen for S0 2 mean 
concentration, which dropped from an annual average of 11.1 
ppb in 1973 to 1.1 ppb in 1991, the last complete year of fol¬ 
low-up. Declines in other pollutant indices were much less 
marked, and most of the cohort experienced many exceed¬ 
ances for the exceedance frequency indices, even in the later 
years of the study. For example, in 1991, two-thirds of the co¬ 
hort experienced one or more months when PM 10 exceeded 
100 pg/m 1 . 

To determine if the elevated risks of CRC mortality as as¬ 
sociated with PMjoQQQ) might be due mainly to levels experi¬ 
enced early in the time period, the proportional hazards re¬ 
gression was rerun using a moving average of PM 10 (100) 
averaged for 4 yr prior to the r isk set. The RR for an IQR in¬ 
crease decreased slightly: RR = 1.15 (95% Cl: 1.00, 1.32). 

Time Spent Indoors 

We reran refined models using adjusted outdoor ambient con¬ 
centrations obtained by applying an indoor/outdoor adjust¬ 
ment factor to mean estimated ambient concentrations ac¬ 
cording to time spent indoors as reported by season for each 
study participant in 1977. Indoor adjustment factors of 0.7 for 
PM10. 0.5 for 0 3 , and 0.6 for SO z were used (17). For N0 2 , In- 
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door sources were not measured until 1987; hence, adjusted 
N0 2 could not be used in mortality analyses (10). The results 
of using adjusted ambient concentration were consistent with 
those reported for unadjusted mean concentrations of PM 10 
with CRC mortality in both sexes, 0 3 with lung cancer mortal¬ 
ity in males, and S0 2 with lung cancer mortality in males and 
in females. 

Restriction of Distance from Monitoring Station 

All interpolations were restricted to be within 31.25 miles of a 
monitoring station. We conducted sensitivity analyses by fur¬ 
ther restricting subjects to those who lived at least 80% of 
their months within 10 miles of a monitoring station for PM lD 
or within 20 miles of a monitoring station for 0 3 . These were 
distances deemed by the Environmental Protection Agency 
(EPA) to be acceptable quality for interpolation (6), Approxi¬ 
mately 63% of subjects met these criteria for PMm and 87% 
for O, When analyses were restricted to these subgroups, the 
RRs relating CRC mortality in both sexes with PM m and iung 
cancer mortality in males with 0 3 were within 6% of previous 
values. The RR of lung cancer mortality with PM 10 was 20% 
smaller for males. 

Comparison of Current with Previous 
Findings from This Cohort 

Previous published findings from analyses on this cohort showed 
no associations between all natural cause mortality and air pol¬ 
lutants during the First 10 yr of follow-up, whereas present anal¬ 
yses using 15 yr of follow-up do show associations (3). Com¬ 
bined male and female RRs of 1.0 or slightly less were observed 
previously The particulate pollutant used previously was days 
per year in excess of 200 pg/m 3 of TSP (TSP[200]). We used 
TSP(20G) averaged over the baseline period 1973-1977 and com¬ 
pared the sex-specific regression coefficients of TSP(200) as asso¬ 
ciated with all natural cause mortality for 10 yr of follow-up ver¬ 
sus 15 yr (1977-1992). The regression coefficients for the two 
time periods were in close agreement. Similar results were ob¬ 
tained for PM| 0 (100). When the model used previously was 
rerun and sex-specific RRs were calculated for an IQR of 
TSP(200) of 33 d/yr, we observed RRs of 1.13 (95% CT 0.97, 
1.31) for males and 0.91 (95% Cl: 0.80, 1.04) for females. Prom 
Table 3, we see that current sex-specific results for PM, 0 were 
very close to these. Sex-specific results were not presented pre¬ 
viously for all natural cause mortality since we had not at that 
time noticed the marked differences between sexes in exposure 
due to males spending considerably more time outdoors than 
females. Further discussion of discrepancies between sexes is 
given below. 

DISCUSSION 

Summary 

In this report, we have found long-term ambient concentra 
tions of PMio tD be associated with increased risks of all natu¬ 
ral cause mortality in males, mortality with any mention of 
nonmalignant respiratory causes in both sexes, and lung can 
cer mortality in males. Long-term ambient concentrations of 
0 3 were even more strongly associated with increased risk of 
lung cancer mortality in males, and SO z was independently as¬ 
sociated with increased risk of lung cancer mortality in both 
males and females. 

In reporting observed associations between specific causes 
of death and individual monitored pollutants, we realize that 
differences could be due to measurement error and that some 
of the observed effects could be the result of correlations 


among pollutants. Also, we were limited to pollutants that had 
been widely monitored over long time periods and recognize 
that the true toxic agent may not have been monitored. 

New Results versus Six Cities and ACS Studies 
Comparisons of current results with published findings from 
the Six Cities and the ACS Studies are somewhat limited be¬ 
cause updated estimates of fine particles less than 2.5 p,m in 
diameter (PM 2 , 5 ) are not yet available on our cohort. The 
ACS Study has published results only pertaining to PM 25 and 
SO, (2). Although the Six Cities Study addressed all of the 
pollutants that We have addressed, theii ability to differentiate 
between specific pollutants was limited due to only having six 
different air pollution sites, those of the six cities. Also, there 
may have been an inadequate range of 0 3 across these cities 
to show effects of 0 3 (1). We compare our results by cause of 
death below with the other two long-term prospective studies 
that have measured respirable particles. For ease of compari¬ 
son with our results given in Table 3, we have converted RRs 
cited by the other studies to correspond to IQR Increases used 
in our study. For the sake of abbreviation, it will be under¬ 
stood that all the RRs cited below are for the IQR increases 
given in Table 3. Results from the Six Cities Study were not 
given sex-specific by pollutant. Comparisons with older stud¬ 
ies that did not measure respirable particles or with cross-sec¬ 
tional studies will not be made. These studies vvere thoroughly 
reviewed in the EPA particulate criteria document (18). Most 
studies show positive associations between long-term ambient 
levels of particulate pollutants and mortality. 

All Natural Cause Mortality 

The Six Cities Study showed associations between all natural 
cause mortality and PM l0 , PM 2 5 , and SO,. The RRs for mean 
concentrations of PM 10 and SO, were 1.23 (95% Cl: 1.07, 1.40) 
and 1.09 (95% CL 1 03, 1 15), respectively. The ACS Study 
showed associations between all natural cause mortality and 
PM 2 5 and SO,. The RR for SO, mean concentration for males 
was 1.02 (95% Cl: 1.01, 1.03) and for females was 1.03 (95% 
Cl: 1.01, 1.04). Our associations of all natural cause mortality 
with PM 10 and SO., for males were positive but smaller than 
for the Six Cities Study and were negative for females. Our 
observed RR for SO, for males was larger than that observed 
in the ACS Study. 

Cardiopulmonary Mortality 

Both the Six Cities and the ACS studies showed positive asso¬ 
ciations between cardiopulmonary deaths and pollutants stud¬ 
ied. The Six Cities Study did not give pollutant-specific RRs 
for Specific causes of death but only the RRs for the most ver¬ 
sus the least polluted city. The RR for the most versus the 
least polluted city was 1.37, somewhat larger than the RR we 
observed for males for PM )0 . The difference in PMi 0 between 
the most and the least polluted cities was 28.3 pg/m 3 (compa¬ 
rable to our IQR): that for S0 4 was 7.5 p.g/m 1 (considerably 
more than our IQR). The ACS Study gave an RR for S0 4 for 
males of 1.03 (95% Cl: 1.01, 1.04) and for females of 1.05 
(95% Cl: 1.03, 1.07). Our data showed a positive association 
for SO, comparable to the ACS Study for males but negative 
associations for PM, 0 and SO, for females (Table 3) 

Any Mention of Nonmalignant Respiratory (CRC) Mortality 

Any mention of nonmalignant respiratory disease as an un¬ 
derlying or contributing cause of death was not studied in the 
other two studies. Our data showed only PMm to have a 
strong association with this outcome for both sexes that re- 
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mained strongly positive when other pollutants were added 
one at a time in two-pollutant models. 

Nonmalignant Respiratory Disease as 
Underlying Cause of Death 

This cause of death as an underlying (not contributing) cause 
of death was considered in the Six Cities Study, but not the 
AC5 Study. The Six Cities Study found no association with 
any pollutant (1). We found a positive, though nonsignificant, 
association with PM l0 in males with an RR of 1.15 {95% Cl: 
0,81, 163). This nonsignificant result may be due to inade¬ 
quate power since there was a small effect and small number 
of cases for this outcome. 

Lung Cancer Mortality 

The Six Cities Study showed an RR for lung cancer mortality 
for the most polluted versus the least polluted city of 1.37 
(95% Cl: 0.81, 2.31) For our IQR of S0 4 , the ACS Study gave 
an RR for males of 1.43 (95% Cl: 1.13. 1.81) and for females 
of 1.17 (95% Cl: 0.80, 1.72). The ACS Study also showed posi¬ 
tive associations between lung cancer mortality and Pl\T 25 . 
However, these associations became negative in never smok¬ 
ers for both sexes. Our data showed a positive association be¬ 
tween PM[,) and lung cancer mortality for females, and our 
strongest association of PMio was with lung cancer mortality 
in males, which remained positive and strong in never smok¬ 
ers. This association also remained positive and strong in two- 
pollutant models with other pollutants. Because of inadequate 
lag time, SO. ; has not been studied in our data with respect to 
this outcome, SO^ data being only available since 1977 in Cali¬ 
fornia. However, SO z was significantly associated with increased 
risk of lung cancer mortality in both sexes, and the effects re¬ 
mained stable in two-pollutant models with other pollutants. 
Long-term averages of SO E and S0 4 were fairly highly correlated 
in our cohort. Ambient levels of 0 3 above 100 ppb were more 
strongly associated with lung cancer mortality in males than were 
any of the other pollutants, but this was not true for females. The 
Six Cities Study lacked the high levels of 0 3 experienced in 
Southern California; the ACS Study did not study O s . 

Differences between Sexes 

Findings were consistent between sexes for the other two 
studies. This was only true in our cohort for associations be¬ 
tween PMia and CRC mortality and between SO2 and lung 
cancer mortality Notable differences between males and fe¬ 
males were observed for lung cancer mortality as associated 
with PM 10 and 0 3 , where males showed much stronger associ¬ 
ations than females. Males In our study spent significantly 
more time outdoors than females, and stronger associations of 
PM 1(t with CRC mortality were observed for males who spent 
more time oucdoors. Thus, differences between males and fe¬ 
males in exposure to ambient pollutants may partially exptaln 
these discrepancies (14). We also noted several differences be¬ 
tween males and females in exposures to other pollutants, such 
as tobacco smoke and occupational sources (Table 1). There 
may be other reasons for the observed differences, and this 
needs to be explored further in future research with a larger 
number of cases. Future research on this cohort pertaining to 
PM 2 5 may find fewer discrepancies between sexes since in¬ 
door and outdoor levels of PM? s tend to correlate highly (2), 
thus diminishing differences between sexes in exposure due to 
males spending more time outdoors. 

Tobacco Smoke Exposure 

Though the effects of tobacco smoke exposure were signifi¬ 
cant for males, the effects were small in this cohort. This may¬ 


be because no reported current smokers were included in 
analyses and the average length of time since past smokers 
had last smoked, estimated from age of joining the church, 
was 27 yr. Only 7.3% of past smokers could be classified as 
possibly smoking within the past 5 yr. 

Results of this study suggest that the associations between 
air pollutants and mortality (both CRC and lung cancer) are 
not due to inadequate control of tobacco smoke exposure be¬ 
cause of the following reasons: ( 1 ) the associations persisted in 
stratified analysis of past smokers and never smokers and re¬ 
mained similar when pack-years of past smoking were con¬ 
trolled for as well as when dummy variables were used to rep¬ 
resent past smoking and estimated years since last smoked; 
and (2) the associations persisted in models that controlled for 
years worked with a smoker and years lived with a smoker. It 
is unlikely that underreporting of current smoking could ac¬ 
count for the observed associations, since our sensitivity anal¬ 
yses indicated that if underreporting were nondifferential, our 
reported RRs would underestimate the true RR, Even if un¬ 
derreporting were differential, reported RRs would overesti¬ 
mate the true RR by not more than 15%. We have no reasons 
to suspect differential underreporting in this cohort. 

Respiratory Disease and Mortality 

The positive association we have observed between death 
with any mention of nonmalignant respiratory disease and 
PM 10 is consistent with earlier findings on this cohort, which 
observed increased prevalence of asthma, chronic bronchitis, 
or emphysema in 1977 with long-term ambient concentrations 
of TSP (19) In current analyses, prevalence of these respira¬ 
tory symptoms in 1977 was found to be a very strong predictor 
of CRC mortality. Recently published analyses pertaining to 
PMjq and respiratory disease on this cohort showed PM 10 to 
be related to development nf overall respiratory symptoms 
and chronic productive cough (4, 20). Elevated long-term con¬ 
centrations of PM 10 were also found to be associated with dec¬ 
rements in lung functions for males, especially those whose 
parents had asthma, chronic bronchitis, emphysema, or hay fe¬ 
ver, and with increased peak expiratory flow lability in both 
sexes (21), Particulate pollutants have also been found to be 
associated with asthma, chronic bronchitis, and emphysema as 
well as decrements in lung function in cross-sectional data 
from U.S. cities (22-24). 

Measurement Error 

We have unknown amounts of measurement error in both the 
estimated long-term ambient concentrations of pollutants as 
well as other covariates. This could bias our estimates of pol¬ 
lutant regression coefficients. One source of measurement er¬ 
ror is that due to interpolating ambient concentrations from 
fixed site monitoring stations to zip code centroids of work 
and home locations of study participants. When cumulated 
over a 2-vr period, correlations between interpolated and ac¬ 
tual monitored values of mean concentration were high—0.88 
for PM 10 , 0.87 for 0 :{ , 0.84 for S0 4 , 0.64 for SO a , and 0.92 for 
NOj (4, 7 25, 26). When analyses were restricted to subjects 
living at least 80% of their months within distances of a moni¬ 
toring station deemed by the EPA to have an acceptable qual¬ 
ity of interpolation, results consistent with those reported above 
without restriction were obtained. Another source of mea 
surement error is our use of indirect estimates of PM| 0 prior to 
1987 when monitored PM| 0 first became widely available 
throughout the state of California. During a 2-yr period, both 
PM|„ and TSP were corvsistantly monitored at a number of sta¬ 
tions throughout California, and these data were used to iden¬ 
tify site- and season-specific regression equations from which 
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PM ln concentrations could be predicted from TSP concentra¬ 
tions (4). For this 2-yr period, we evaluated the precision of 
our interpolation algorithms for predicting PM 10 concentra¬ 
tions at each monitoring site using the PM 1C , measurements 
from the surrounding monitoring sites and also using the PM J0 
values estimated From the TSP measurements from the same 
surrounding monitoring sites. The correlation between mea¬ 
sured mean PM 10 and interpolated mean PM ]0 was 0.88 and 
0.86, respectively, using measured PM K! and PM ]0 estimated 
from TSP values from the surrounding monitoring sites. The 
correlation between measured PM 10 (100) and interpolated 
PMjo(IOO) was 0.79 and 0.67, respectively, using measured 
PM 1B and PM 1C , estimated from TSP values from the surround¬ 
ing sites. This indicates that use of TSP to predict mean PM, p 
resulted in minimal Idss of precision compared with use of 
measured PM 10 Use of TSP as a predictor for PM ]n (100), 
however, for the years prior to 1987 resulted in some addi¬ 
tional measurement error. The consistency in our findings re¬ 
garding health effects of the different indices of PM W suggests 
that this source of measurement error may not be seriously af¬ 
fecting our results. 

Use of Ambient Concentration Estimates for Pollutants 

An Important limitation of assessing associations between 
long-term ambient concentrations of air pollutants and mor¬ 
tality is that our estimates are of outdoor ambient concentra¬ 
tions, not personal exposure. We have dealt with the lack of 
personal exposure measurements in two ways: (/) by offering 
self-reports of personal exposure to other sources of particles, 
such as past smoking, passive smoking, occupational, etc., as 
covariate candidates for the health effects models; and (E) by 
adjusting mean ambient concentrations of all pollutants ac¬ 
cording to time spent indoors. The results of using adjusted 
ambient mean concentrations were consistent with those re¬ 
ported for unadjusted mean concentrations. 

Conclusions and Future Research Needed 

Our results indicate that long periods of residence and work 
location in areas of high ambient air pollution are associated 
with increased mortality. We observed associations of PM 10 
with all natural cause mortality in males and with nonmalig- 
nant respiratory mortality in both sexes. Lung cancer mortal¬ 
ity was associated with PM ln and 0 3 in males and with S0 2 in 
both males and females. These results do not appear to be due 
to confounding by any of a large number of measured risk fac¬ 
tors In this cohort of nonsmoking California adults. 

Although findings of the relationship between PM 10 and 
mortality from our study are qualitatively in agreement with 
those from the Six Cities and ACS Studies, the discrepancies 
among these studies reported above, with respect to specific 
particulate pollutants and specific causes of mortality, need to 
be further explored (1, 2). This would require updating PM i 5 
estimates for our cohort and adding other pollutants to the 
ACS Study, as well as studying contributing nonmalignant re¬ 
spiratory causes of death for the other two studies. 
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APPENDIX 

Additional Factors Evaluated for Confounding or 
Effect Modification 

The factors evaluated for additional models were the follow¬ 
ing; whether father or mother had cancer; whether subject re¬ 
ported having high blood pressure at baseline; a physical exer¬ 
cise scale composed from a number of questions on leisure 
time and work physical activity and categorized as rone, light. 


moderate, or heavy; whether subject currently used alcoholic 
beverages at baseline; times per week and dosage per week of 
antioxidant vitamins A, C, or E, as consumed in pills or multi¬ 
vitamin tablets; an index of magnesium intake; and times per 
week subject used cooked green vegetables, green salads, to¬ 
matoes, fruit or fruit juice, beans or lentils, and eggs Various 
indices of meat, poultry, and fish consumption were also eval¬ 
uated as well as years of the SDA lifestyle as measured by age 
at baptism into the SDA Church. Covariates that could mod¬ 
ify exposure to ambient air pollutants or exacerbate health ef¬ 
fects of such exposure were evaluated for interaction with the 
pollutant using the cross-product term of each such variable 
with the pollutant. Covariates evaluated for interaction in¬ 
ch] ded hours outdoors, computed as a weighted average of 
summer and winter seasons in 1977, antioxidant food intake as 
measured by the vegetable/fruit index; antioxidant vitamin in¬ 
dex from supplements; fish consumption; and whether subject 
had heart attack, high blood pressure, stroke, or diabetes as 
reported at baseline. None of these showed significant interac¬ 
tions with PM io or effect modifications in stratified analyses. 
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